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ABSTRACT 

We present deep UBVRI CCD photometry for the young open star clusters Tr 1 and 
Be 11. The CCD data for Be 11 is obtained for the first time. The sample consists 
of ^ 1500 stars reaching down to y 21 mag. Analysis of the radial distribution of 
stellar surface density indicates that radius values for Tr 1 and Be 11 are 2.3 and 1.5 
pc respectively. The interstellar extinction across the face of the imaged clusters region 
seems to be non-uniform with a mean value of E{B — V) — 0.60±0.05 and 0.95±0.05 
mag for Tr 1 and Be 11 respectively. A random positional variation of E{B — V) is 
present in both the clusters. In the cluster Be 11, the reason of random positional 
variation may be apparent association of the HII region (S 213). The 2MASS JHK 
data in combination with the optical data in the cluster Be 11 yields E{J — K) = 
0.40±0.20 mag and EiV — K) ~ 2.20±0.20 mag. Colour excess diagrams indicate a 
normal interstellar extinction law in the direction of cluster Be 11. 
The distances of Tr 1 and Be 11 are estimated as 2.6±0.10 and 2.2±0.10 Kpc respec- 
tively, while the theoretical stellar evolutionary isochrones fitted to the bright cluster 
members indicate that the cluster Tr 1 and Be 11 are 40±10 and 110±10 Myr old. The 
mass functions corrected for both field star contamination and data incompleteness 
are derived for both the clusters. The slopes 1.50±0.40 and 1.22±0.24 for Tr 1 and Be 
11 respectively are in agreement with the Salpeter's value. Observed mass segregations 
in both clusters may be due to the result of dynamical evolutions or imprint of star 
formation processes or both. 

Key words: Star clusters - individual: Tr 1 and Be 11 - star: Interstellar extinction, 
luminosity function, mass function, mass segregation - HR diagram. 



1 INTRODUCTION 

Young open star clusters in a galaxy provide valuable infor- 
mation about star formation processes and are key objects 
for the galactic structure and evolution. For such studies, 
a knowledge of cluster's parameters like distance, age, red- 
dening and stellar content is required which can be derived 
from the colour-magnitude (CM) and colour-colour (CC) 
diagrams of star clusters. In addition to this, the distribu- 
tion of stellar masses at the time of cluster formation is of 
fundamental importance to analysis related to evolution of 
galaxies. The initial mass function (IMF) also plays an im- 
portant role in understanding the early dynamical evolution 
of star clusters, because it is a fossil record of the very com- 
plex process of star formation and provides an important 
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link between the easily observable population of luminous 
stars in a stellar system and the fainter, but dynamically 
more important, low mass stars. Another related problem is 
the mass segregation in star clusters in which massive stars 
are more concentrated towards the cluster center compared 
to low mass stars. It is not clear whether the mass segrega- 
tion observed in several open clusters is due to dynamical 
evolution or an imprint of star formation processes itself (cf. 
Sagar et al. 1988; Sagar 2001 and references therein). Thus 
one can say that the young open star clusters are the labo- 
ratories in a galaxy for providing answers to many current 
questions of astrophysics. 

In the light of above discussions, we performed multi- 
colour deep CCD stellar photometry in two young open star 
clusters namely Trumpler 1 (Tr 1) and Berkeley 11 (Be 11). 
The CCD UBVRI observations of Be 11 are presented for 
the first time. The relevant prior informations (taken from 



© 2002 RAS 



2 



R. K. S. Yadav and Ram Sagar 

Table 1. General information about the clusters under study, taken from Mermilliod (1995) 



Cluster 


lAU 


OCL 


1 


b 


Trumpler 


Radius 


Distance 


E{B - V) 


log(age) 








(deg) 


(deg) 


class 


(arcmin) 


(Kpc) 


(mag) 


(ys) 


Trumpler 1 


C0132+610 


328 


128.22 


-1.14 


II 2p 


1.5 


2.6 


0.58 


7.5 


Berkeley 11 


C0417+448 


404 


157.08 


-3.65 


II 2m 


2.5 


2.2 


0.95 


7.7 



Mermilliod (1995)) of these clusters are given in Table 1. 
These clusters are relatively compact objects with angular 
radii less than 3'. Previous studies, observations and data re- 
ductions are described in the next sections. The interstellar 
extinction, other photometric results, luminosity function, 
mass function and mass segregation arc described in the 
subsequent sections. Making use of JHK data with optical 
data, extinction law has also been studied in Be 11. 



2 PREVIOUS STUDIES 

Trumpler 1: It is an extremely concentrated galactic open 
star cluster in Cassiopeia. It lies at the outer edge of the 
Perseus spiral arm. Oja (1966) carried out the proper motion 
study. UBV photoelectric photometry for 43 bright stars 
was presented by Joshi & Sagar (1977) while Phelps & Janes 
(1994) published UBV CCD photometry (l^~18 mag). Mc- 
Cuskey & Houk (1964) studied this cluster photographically 
in UBV system while Steppe (1974) has presented three 
colour RGU photographic photometry. All these studies in- 
dicate that reddening across the cluster is uniform with 
E{B — V) = 0.61 mag, distance estimate is 2630 pc and 
age seems to be ~ 27 Myr. 

Berkeley 11: It is a distant neglected compact young open 
cluster, apparently associated with the faint HII region S 
213. Only UBV photoelectric photometry of 24 bright stars 
has been done by Jackson et al. (1980). On the basis of 
this study, they found that this cluster has members earliest 
photometric type ^ b4 and is thus an extreme Population I 
object. They also found reddening E{B - V) = 0.95±0.06 
mag, distance d = 2.2± 0.2 Kpc and age of the cluster as 
3x10'' yr. 



3 OBSERVATIONAL DATA 

The optical and near-IR JHK data used in the present 
study are described in the following subsections. 

3.1 Optical observations and data reductions 

The UBV Johnson and RI Cousins observations were ob- 
tained, using 2Kx2K CCD system at the f/13 Cassegrain 
focus of the Sampurnanand 104-cm telescope of the State 
Observatory, Nainital during November 2000. Details of the 
observations are given in Table 2. Each pixel of 2048x2048 
size CCD corresponds to a square of size 0".36 on the sky. lu 
order to improve the S/N ratio, the observations were taken 
in binning mode of 2x2 pixel. The entire chip covers a field 
of 12'.3xl2'.3 though because of the smaller filter sizes, we 



could image only ~ 8'. 6x8'. 6 region. The read-out noise for 
the system is 5.3 e~ with a gain of 10 e~/ADU. In each 
passband, only one short but 2 to 3 deep exposures were 
taken for the cluster region so that accurate photometric 
measurements can be obtained for faint stars. For remov- 
ing field star contamination, we also observed field regions 
in the UBV{RI)c passband situated ~ 15' south for both 
the clusters. Fig. 1 show the identification maps of imaged 
cluster and field regions of Tr 1 and Be 11. For calibra- 
tion, we observed 11 Landolt (1992) standard stars covering 
a range in brightness (11.0<y <15.0) as well as in colour 
(0.16<(y - /)<2.08). Flat field exposures ranging from 20 
to 60 sec in each filter were made on the twilight sky. A 
number of biases were also taken during the observing runs. 

The CCD data frames were reduced using computing 
facilities available at the State Observatory, Nainital. Ini- 
tial processing of the data frames were done in the usual 
maimer using the IRAF data reduction package. Different 
cleaned frames of the same field in the same filter were co- 
added. Photometry of co-added frames was carried out us- 
ing DAOPHOT software (Stetson 1987). PSF was obtained 
for each frame using several uncontaminated stars. In those 
cases where brighter stars are saturated on deep exposure 
frames, their magnitudes have been taken only from the 
short exposure frames. Wherever more than one measure- 
ment is available in a passband for a star, the final magni- 
tude is an average of the individual measurements and its 
error is the ALLSTAR error of the average. When only one 
measurement is available, the error is taken to be the output 
of ALLSTAR. 

The photometric calibration equations are determined 

by fitting least square linear regression to the standard 
UBVRI photometric indices as function of observed instru- 
mental magnitudes normalized for 1 second exposure time. 
The following colour equations are obtained for the system. 

A{U -B) = (0.973±0.020)A('u - b)o 
A{B -V) ^ (1.117±0.014)A(6 - v)o 
AiV - R) = (0.886±0.011)A('(; - r)o 
A(V' - /) = (0.992±0.007)A(v - i)o 
AY = Avo - (0.060±0.024)(y - 7) 

where A denotes differential values; (U — B), (B — V), 
{V — R), {V — I) and V are standard values taken from 
Landolt (1992) and {u — b)o, {b~-v)o, {v — r)o {v~i)o and iiq 
are the instrumental CCD aperture magnitudes and colours 
corrected for atmospheric extinction. The atmospheric ex- 
tinction coefficients are 0.59±0.07, 0.37±0.02, 0.29±0.06, 
0.16±0.03 and 0.13±0.02 for U, B, V, R and I respectively. 
The errors in the colour coefficients are obtained from the 
deviation of data points from the linear relation. These equa- 
tions are used to standardize the CCD instrumental mag- 
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Figure 1. Identification maps for tlie cluster and field regions of Tr 1 and Be 11. The (X, Y) coordinates are in pixel units corresponding 
to 0".72 on the sky. East is up and north is left. Filled circles of different sizes represent brightness of the stars. Smallest size denotes 
stars of mag. Open circles in the clusters region represent the cluster size. 



nitudes of both cluster and field regions. To establisii tiie 
local standards, we selected several isolated stars in the ob- 
served regions and used the DAOGROW programme for the 
construction of an aperture growth curve required for deter- 
mining the difference between aperture and profile - fitting 
magnitudes. These differences and differences in exposure 
times and atmospheric extinctions are used in evaluating 
zero points for the reference frames. The zero points are 
uncertain by 0.02 mag in V; ~ 0.01 mag in {B — V), 
{V ~R), {V - I) and ~ 0.03 mag in (U - B). Other factors 
contributing to the photometric uncertainty are described 
recently by Moitinho (2001). Amongst them, the internal 
errors estimated on the S/N ratio of the stars as output of 
the ALLSTAR mainly produce the scatter in the various CC 
and CM diagrams of the clusters. They are given in Table 3 
as a function of brightness for the cluster region. The errors 
become large (>0.1 mag) for stars fainter than V—20 mag, 
so the measurements should be considered unreliable below 



this magnitude. 



The (X, Y) pixel coordinates as well as the V, {U — B), 
{B — V), {V — R) and (V — 7) magnitudes of a sam- 
ple of stars observed in the cluster and field regions of 
Tr 1 and Be 11 are listed in Tables 4 and 5 respectively. 
Only formats of the tables are presented here; the full ta- 
bles are available in the electronic version of the article on 
Synergy, on the open star cluster data base web site at 
http://obswww.unige.ch/wehda/, and also with the authors. 
In order to avoid introducing a new numbering system in Tr 
1, we adopt the numbers from the data base given by Phelps 
& Janes (1994). Stars not observed earlier have a number 
starting with 1001 in Tr 1. 
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Table 2. Log of CCD observations. N denotes the number of stars measured in different passbands. 



Region 


Filter 


Exposure Time 
(in seconds) 


Date 


N 


Trumpler 1 - cluster 


U 


1800x2, 300x1 


19/20 Nov 2000 


524 


f^nnnr, — ^^ ^ '^^'^ df\^ 
L*200U — 




l-^\J\J AO, ^'±U /V J- 




1283 


52000 = +61<*17'20" 


V 


900x3, 120x1 




1451 




R 


480x3, 60x1 




1630 




I 


240x3, 60x1 




1710 


Trumpler 1 - field 


U 


1200x1 


20/21 Nov 2000 


350 


02000 = 0l''35'"40= 


B 


900x1 




915 


52000 = +61'^02'20" 


V 


600x1 




1190 




R 


300x1 




1346 




I 


300x1 




1585 


Berkeley 11 - cluster 


U 


1800x2, 300x1 


18/19 Nov 2000 


295 


aaooo = 04''20'"36'' 


B 


1200x3, 240x1 


5 5 


650 


52000 = +44'^55'58" 


V 


900x3, 180x1 


1 1 


750 




R 


600x3, 120x1 


1 1 


930 




I 


300x3, 60x1 


1 1 


1090 


Berkeley 11 - field 


U 


900x1 


20/21 Nov 2000 


90 


"2000 = 04''20'"34'' 


B 


600x1 


?? 


299 


52000 = +44'*40'08" 


V 


600x1 


)) 


474 




R 


300x1 


)) 


676 




I 


300x1 


II 


723 



Table 3. Internal photometric errors as a function of brightness. 
a is the standard deviation per observation in magnitude. 
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Figure 2. A comparison of present photometry with CCD and 
photoelectric data of Tr 1 and Be 11. For Tr 1, filled circles rep- 
resent the CCD data of Phelps & Janes (1994) and triangles rep- 
resent the photoelectric data of Joshi & Sagar (1977). For Be 
11 filled circles represent the photoelectric data of Jackson et al. 
(1980). 



Magnitude range 


0"!7 


OB 


ay 






<12.0 


0.012 


0.001 


0.005 


0.003 


0.011 


12.0 - 13.0 


0.013 


0.004 


0.008 


0.008 


0.013 


13.0 - 14.0 


0.019 


0.006 


0.010 


0.009 


0.019 


14.0 - 15.0 


0.021 


0.007 


0.014 


0.012 


0.014 


15.0 - 16.0 


0.023 


0.012 


0.015 


0.013 


0.014 


16.0 - 17.0 


0.025 


0.012 


0.017 


0.015 


0.020 


17.0 - 18.0 


0.043 


0.024 


0.026 


0.025 


0.030 


18.0 - 19.0 


0.083 


0.054 


0.055 


0.051 


0.060 


19.0 - 20.0 


0.086 


0.095 


0.102 


0.091 


0.113 


20.0 - 21.0 


0.158 


0.281 


0.356 


0.346 


0.264 



present CCD data for Tr 1 show a constant zero point 
offset of ~ 0.08 mag in AV with that of given by Phelps 
& Janes (1994). A weak linear dependence of A(B — V) is 
observed on brightness. The A(L'' — B) values, on the other 
hand, show decreasing trend upto V ~ 16.0 mag but increase 
for fainter stars. As UBV photoelectric data is in good 
agreement with the present CCD data for both the clus- 
ters, we suspect calibration problem with Phelps & Janes 
(1994) CCD data. 



3.2 Comparison with the previous photometries 

We compare the present data with the published CCD and 
photoelectric data. Table 6 represents the average differences 
in the sense present minus others along with their standard 
deviations. The differences A in mag and colours are plotted 
in Fig 2. Fig 2 and Table 6 indicates that 



3.3 Near - IR data 

The near-IR JHK data are taken from the digital Two 
Micron All Sky Survey (2MASS) available at web site 
http://www.ipac.caltech.edu/2MASS/. 2MASS is uniformly 
scanning the entire sky in three near-IR bands J (1.25 
^m), 77(1.65 ^m) and 7^8(2.17 ^m) with two highly - 
automated 1.3-m telescopes equipped with a three chan- 
nel camera, each channel consisting of a 256x256 array 
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Table 4. CCD relative (X, Y) positions and V, (U - B), (B - V), (V - R) and {V - I) photometric 
magnitudes of few stars, as a sample measured In the cluster and field region of the cluster Tr 1. In the 
cluster region, stars observed earlier have numbering system of Phelps &; Janes (1994) taken from the cluster 
data base, while numbering of stars observed for the first time by us start with 1001 in the column 1. The 

last column represent the photometric membership informations where m and nm represent the member and 
non-member stars. In the field region, stars are numbered in the increasing order of X value. 



Star 


X 


Y 


V 


(U-B) 


(B-V) 


(V-R) 


(V-I) 


Mem 




(pixel) 


(pixel) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 












Tr 1 cluster region 










1001 


97.64 


229.68 


19.80 


* 


1.21 


0.72 


1.49 


nm 


1002 


356.27 


238.04 


18.67 


0.42 


1.16 


0.66 


1.35 


m 


1003 


338.42 


237.33 


19.54 


* 


1.31 


0.78 


1.54 


m 


1004 


309.76 


235.40 


18.04 


0.48 


1.17 


0.66 


1.31 


m 


1005 


155.45 


227.56 


17.66 


0.61 

Tr 1 field region 


1.08 


0.65 


1.26 


nm 


1 


3.32 


151.93 


19.36 


* 


1.20 


0.68 


1.46 




2 


6.93 


238.09 


18.45 


* 


1.23 


0.67 


1.45 




3 


14.91 


479.27 


19.46 


* 


1.21 


0.73 


1.60 




4 


15.22 


219.91 


18.94 


* 


1.37 


0.67 


1.52 




5 


16.93 


360.70 


19.06 


* 


1.30 


0.65 


1.46 





Table 5. CCD relative (X, Y) positions and V, (U-B), (B-V), (V-R) and (V-/) photometric magnitudes 
of few stars, as a sample measured in the cluster and field regions of the cluster Be 11. Stars are numbered 
in increasing order of X value. The last column indicates the photometric member and non-member stars of 
the cluster. 



Star 


X 


Y 


V 


(U-B) 


(B-V) 


(V-R) 


(V-I) 


Mem 




(pixel) 


(pixel) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 










Be 11 cluster region 












142 


311.95 


318.82 


14.78 


0.61 


0.99 


0.48 


1.04 


m 


143 


312.78 


178.07 


16.97 


0.77 


1.30 


0.63 


1.33 


m 


144 


313.98 


339.75 


15.11 


0.43 


0.91 


0.45 


0.97 


m 


145 


316.18 


128.17 


18.53 


* 


1.54 


0.77 


1.65 


nm 


146 


317.11 


151.06 


19.43 
Be 11 field region 


* 


1.81 


0.92 


1.88 


m 


1 


14.00 


268.90 


18.02 


* 


2.10 


1.15 


2.40 




2 


29.36 


577.23 


17.64 


0.60 


1.10 


0.62 


1.29 




3 


30.76 


610.82 


20.98 


* 


* 


0.84 


1.65 




4 


36.19 


558.98 


20.62 


* 


* 


0.77 


1.56 




5 


37.84 


701.05 


14.40 


0.63 


1.10 


0.50 


1.06 





of HgcdTe detectors. The photometric uncertainty of the 
data is < 0.155 mag with Ks ~ 16.5 mag photomet- 
ric completeness. The Ks magnitudes are converted into 
K magnitude following Persson ct al. (1998). Further do- 
tails about the 2MASS and 2MASS data are available 
at http://www.ipac.caltech.edu/2mass/releases/second/doc/ 
exphup .html. The JHKs data arc available for 200 stars in 
the observed area of Be 11 cluster. Among them only 179 
stars are common with our optical data. 



4 DATA ANALYSIS 

4.1 Cluster radius and radial stellar surface 
density 

The first step to determine the cluster radius is to find its 
center in the image. The center of the cluster is determined 
iteratively by calculating average X and Y positions of the 
stars within 300 pixels from an eye estimated center, until 



they converged to a constant value. An error of a few tens of 
pixels is expected in locating the cluster center. The (X, Y) 
pixel coordinates of the cluster centers obtained in this way 
are (345, 415) and (326, 432) for Tr 1 and Be 11 respectively. 
The corresponding equatorial coordinates are given in Table 
2. For determining the radial surface density of stars p(r) in 
a cluster, the imaged area has been divided into a number of 
concentric circles with respect to the above estimated cen- 
ter, in such a way that each zone contains a statistically 
significant number of stars. The number density of stars, p,, 
in the i"* zone has been calculated as = where Ni is 

the number of stars and Ai is the area of the i*** zone. The 
density versus radius plots for Tr 1 and Be 11 are shown in 
Fig 3. A clear radius-density gradient present in Fig 3 con- 
firms the existence of clustering. Following Kaluzny (1992), 
we describe the p(r) of an open cluster as: 



" l + (r/re)^ ' 
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Table 6. Comparison of our photometry with others for the cluster Tr 1 and Be 11. The difference (A) is 
always in the sense present minus comparison data. The mean and standard deviations in magnitude are 
based on N stars. Few deviated points arc not included in the average determination. 



Cluster Comparison data V range < AV > < A{B - V)> < A{U - B) > 

Mean±(T(N) Mcan±(T{N) Mcan±o-(N) 



ir i r fiGips qz janes i^iyy^j 


< 




n nc;_i_n c\oi'7\ 
U.UOitU.Uzl^ ( ) 


U.UoztU.Uli^o ) 


U. iDitU.UD^D ) 




14.0 


- 15.0 


0.06±0.02{17) 


0.04±0.02(17) 


0.07±0.07(14) 




15.0 


- 16.0 


0.05±0.06(19) 


0.02±0.02(16) 


0.02±0.07(15) 




16.0 


- 16.5 


0.08±0.02(12) 


0.01±0.03(12) 


0.03±0.05(12) 




16.5 


- 17.0 


0.08±0.03(10) 


0.03±0.03(10) 


0.11±0.07(10) 




17.0 


- 17.5 


0.08±0.03{14) 


0.01±0.08(13) 


0.13±0.10(11) 




17.5 


- 18.0 


0.08±0.04{37) 


-0.03 ±0.05(37) 


0.12±0.09(29) 




18.0 


- 18.5 


0.08±0.05(30) 


-0.03 ±0.09(30) 


0.10±0.11(13) 




18.5 


- 19.0 


0.07±0.03(43) 


-0.02 ±0.07(43) 




Joshi & Sagar (1977) 


12.0 


- 13.0 


-0.01 ± 0.14(3) 


0.03±0.04(4) 


0.03±0.07(4) 




13.0 


- 14.0 


0.00±0.06(10) 


0.05±0.04(12) 


0.04±0.07(12) 




14.0 


- 15.0 


0.03±0.09(7) 


0.07±0.08(17) 


0.03±0.08(17) 


Be 11 Jackson et al. (1980) 


13.0 


- 14.0 


-0.05 ± 0.02(6) 


0.00±0.01(6) 


-0.03 ±0.01(6) 




14.0 


- 15.0 


-0.03 ± 0.01(7) 


-0.03 ±0.01(7) 


-0.02 ±0.03(7) 




15.0 


- 16.0 


-0.03 ± 0.02(5) 


-0.02 ±0.03(5) 


-0.02 ±0.04(5) 



100 300 300 

' \ ' \ ' I 

5,0 - 

Tr 1 




100 yoo 

Radius (Pixel) 

Figure 3. Radial density profile for Tr 1 and Be 11. The length 
of the errorbar denotes errors resulting from sampling statistics 
(=— ^ where N is the number of stars used in the density esti- 
mation at that point). Dotted curves represent fitted profile and 
arrows represent the level of field star densities. 



where the cluster core radius rc is the radial distance at 
which the value of p(r) becomes half of the central density 
/o. We fit this function to the observed data points of each 
cluster and use minimization technique to determine 
and other constants. As can be seen in Fig 3, the fitting of 
the function is satisfactory. The values of core radii derived 
in this way are 85±10 and 180±10 pixels for Tr 1 and Be 11 
respectively. Fig 3 also separate the cluster region from the 
surrounding field region. The field star density thus obtained 
are 2.2x10"^ and 1.0x10"^ stars/pixel^ for Tr 1 and Be 11 



respectively. The radial distribution of stars in Tr 1 and Be 
11 indicates that the extent of the cluster is about 250 and 
200 pixels respectively which correspond to ~ 3.0 and 2.4 in 
arcmin. The cluster sizes are thus a few times larger than 
the corresponding core sizes which is in agreement with the 
findings of Nilakshi et al. (2002). 



4.2 Apparent colour-magnitude diagrams of 
cluster and field regions 

The apparent CM diagrams generated from the present data 
for the Tr 1 and Be 11 clusters and their field regions are 
displayed in Fig 4. The CM diagrams extend up to 1/ ~ 21 
mag except in V , {U — B) diagram where it is only up to V 
~ 18 mag. A well defined cluster MS contaminated by field 
stars is clearly visible in all CM diagrams. The field star con- 
tamination increases with decreasing brightness. The cluster 
sequence fainter than V ^ 17 mag have larger scatter. This 
may be due to photometric errors as well as field star con- 
tamination. It is difficult to separate field stars from the 
cluster members only on the basis of their closeness to the 
main populated area of the CM diagrams, because field stars 
at cluster distance and reddening also occupy this area (see 
Fig 4) . For the separation of cluster members from the field 
stars, precise proper motion and/or radial velocity measure- 
ments of these stars are required. In the absence of such 
data, we use photometric criterion for separating obvious 
field stars. A star is considered as a non-member if it lies 
outside the cluster sequence in at least one CC or CM di- 
agrams. From the V , {V — I) diagram of the field region, 
statistically expected number of field stars among the pho- 
tometric cluster members has been given in Table 7. The 
frequency distribution of the field star contamination in dif- 
ferent part of the CM diagram can be estimated from the 
Table 7. It is thus clear that all photometric probable mem- 
bers can not be cluster members and non-members should 
be subtracted in the studies of cluster MF etc. However, 
probable members located within a cluster radius from its 
center can be used to determine the cluster parameters, as 
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Figure 4. The V , (U — B); V, (B — V) and V, {V — I) diagrams for the stars observed by us in Tr 1 and Be 11 cluster regions and the 
V, {V — I) CM diagram of the corresponding field regions. Dotted lines represent the blue and red envelope of the cluster MS. 



they have relatively less field star contamination and this 
has been done in the sections to follow. 



4.3 Interstellar extinction in the direction of 
clusters 

To estimate interstellar extinction in the direction of the 
clusters, we plot in Fig 5 {U — B) versus {B — V) diagrams 
of the sample stars. Adopting the slope of reddening line 
E{U - B)/E{B - V) as 0.72, we fit the intrinsic zero-age 
main-sequence (ZAMS) given by Schmidt-Kaler (1982) to 
the MS stars of spectral type earlier than AO. This gives a 
mean value of E{B — V) = 0.60±0.05 mag for the cluster Tr 
1 and 0.95±0.05 mag for the cluster Be 11. Our reddening 
estimates for the imaged region agree fairly well with the 
values estimated earlier by others (see Table 1). 

For determining the nature of interstellar extinction law 



in the direction of clusters, we used the stars having spec- 
tral type earlier than AO. This has been selected from their 
location in the {U — B) versus [B — V) and apparent CM di- 
agrams which reveals that bright stars {y< 15.0 mag) with 
(B - l/)<0.90 mag in Tr 1 and with {B - V) < 1.30 mag 
in Be 11 are wanted objects. For these stars, their intrinsic 
colours have been determined using either the spectral type 
(available only for 15 stars in Be 11) taken from the open 
cluster data base (cf. Mermilliod 1995) or UBV photomet- 
ric Q-method (cf. Johnson & Morgan 1953; Sagar & Joshi 
1979) and the calibrations given by Caldwell et al. (1993) 
for (U - B)o, {V - R)o and (V - I)o with (B - V)o. The 
mean values of the colour excess ratios derived in this way 
are listed in Table 8 for both the clusters. They indicate 
that the law of interstellar extinction in the direction of the 
clusters under discussion is normal. 
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Table 7. Frequency distribution of the stars in the V, {V — I) diagram of the cluster and field regions. 

S'lid ^-R denote the number of stars in a magnitude bin blueward, along and redward of the 
cluster sequence respectively. Nc (diff'erence between the Ns value of cluster and field regions) denotes the 
statistically expected number of cluster members in the magnitude bin. 



Be 11 

V range Cluster region Field region Cluster region Field region 

Nb Ns Nr Nb Ns Nr Nc Ng Ns Nr Nb Ns Nr Nc 



13 - 14 





8 


8 


1 


3 


1 


5 





10 














10 


14 - 15 





20 


25 





3 


21 


17 





16 


4 





4 





12 


15 - 16 





22 


40 


1 


5 


60 


17 





26 


10 





8 


1 


18 


16 - 17 


3 


37 


53 


4 


25 


80 


12 





48 


14 


6 


9 


10 


39 


17 - 18 


18 


97 


75 


24 


74 


74 


23 





68 


15 


22 


16 


14 


52 


18 - 19 


52 


131 


71 


80 


110 


52 


21 


14 


142 


19 


33 


22 


14 


120 


19 - 20 


130 


139 


51 


105 


110 


28 


29 


45 


159 


14 


68 


34 


14 


125 




V.{3 V) = 0.95 
E(U-B) = 0.6a 




(B-V) 



E(B-V) 



Figure 5. The {U — B) versus (B — V) diagrams of the stars in 
cluster region observed by us for Tr 1 and Be 11. The continuous 
straight line represents slope 0.72 and direction of the reddening 
vector. The dotted curve represents the locus of Schmidt - Kaler's 
(1982) ZAMS fitted for the marked values of colour excesses. 

4-3.1 Presence of non-uniformity in E(B — V) 

In order to see the extent of non-uniform extinction in the 
clusters under study, we plot histograms of E{B — V) in 
Fig 6. They indicate that there is a range of ~ 0.4 mag 
in E{B — V) values of both the cluster stars with peak 
around E{B - V) = 0.65 and 0.95 mag for Tr 1 and Be 
11 respectively. To determine the presence of non-uniform 
extinction in a cluster, we calculate the value of AE{B — V) 
= E{B - V)max-E{B ~ V)min, where E{B ~ V)max and 
E{B — V)min are determined on the basis of respectively, the 
five highest and five lowest E{B—V) values of the MS cluster 
members. In this way we can find the values of AE{B — V) 
and these values are listed in Table 9. As the factors other 
than non-uniform extinction like stellar evolution, stellar du- 
plicity, stellar rotation, difference in chemical composition. 



Figure 6. Solid and dotted lines represent the histograms of 
E(B — V) for Tr 1 and Be 11 respectively. They indicate presence 
of non - uniform extinction in both clusters. 

dispersion in ages, distances, and uncertainty in photometric 
data (cf. Burki 1975, Sagar 1987, Yadav & Sagar 2001) can 
produce a maximum dispersion in E{B — of ~ 0.11 mag 
for MS members, we consider the presence of non-uniform 
extinction in the both cluster regions, as the observed val- 
ues of AE{B — V) for their MS stars are much more than 
0.11 mag. In order to see whether it is due to presence of 
varying amount of matter inside them or due to any other 
regions, we studied below the variation of E{B — V) with 
spatial position of stars in the clusters. 

4-3.2 Spatial variation of E{B ~ V) 

To study the spatial variation of reddening in terms of E{B — 
V) across the cluster region, we divide the cluster field into 
equal area of small boxes of size I'.Oxl'.O. The positional 
variation of E{B — V) is shown in Table 10 for both the 



© 2002 RAS, MNRAS 000, 1-15 



A deep UBVRI CCD photometric study of the open clusters Tr 1 and Be 11 9 



Table 8. A comparison of the colour excess ratios with E(B — V) 
for star clusters with the corresponding values for the normal 
interstellar extinction law given by Cardelli et al. (1989). 



Objects 


E{U-B) 
E(B-V) 


E(y~R) 

E(B~V) 


E(V-I) 
E(B-V) 


Normal interstellar 


0.72 


0.60 


1.25 


Tr 1 


0.72±0.03 


0.54±0.02 


1.13±0.05 


Be 11 


0.71±0.02 


0.45±0.02 


1.04±0.04 



Table 9. The values olE{B-V)miu, iB-V)max a.ndAE(B-V) 
(see text). 



Cluster 


E{B — V)min 


E{B — V)max 


/^E{B - V) 




(mag) 


(mag) 


(mag) 


Tr 1 


0.55 


0.85 


0.30 


Be 11 


0.91 


1.06 


0.15 



clusters. An inspection of Table 10 indicates that E{B — 
V) does not show any systematic variation with position in 
both the clusters. However, it varies randomly which may 
be due to the random distribution of gas and dust within 
the clusters. 
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4-3.3 Interstellar extinction in near - IR 

In order to derive the interstellar extinction for Be 11 in 
near - IR region, we have combined optical and infrared 
data. There are 47 common stars within the cluster radius. 
We construct in Fig. 7 (J — K) vs {V — K) diagram for the 
stars within the cluster radius and fit a ZAMS for metallic- 
ity Z = 0.02 taken from Schaller et al. (1992). This yields 
E{J-K) = 0.40±0.20 mag and E{V-K) = 2.20±0.20 mag, 
which corresponds to a ratio — 0.18±0.30 in good 

agreement with the normal interstellar extinction value 0.19 



1.6 - 

K(J-K) ^ 0.40 + 0,30 

E(v-K) = g.aoio.ao 



1.2 - 



0.8 - 



0.4 - 




2 3 4 5 6 

V-K 



Figure 7. The (J — K) versus {V — K) colour - colour diagram 
of all the stars which are common in V and JHK data within 
the cluster radius for the cluster Be 11. The solid line is a ZAMS 
fitted for the marked values of colour excesses. 



Figure 8. The plot of E{U - B), E{B - V), E(V - R) E(V - 1), 
E{V - H) and (V - K) against E(V - J) for Be 11. Solid line in 
each diagram represents least square linear fit to the data points. 
The values of correlation coefficients are shown in the diagram. 

suggested by Cardelli et al. (1989). However, scattering is 
larger due to the error size in JHK data. 

4-3.4 Extinction Law in Be 11 

In order to see the nature of extinction law in Be 11, we plot 
the colour excess E{U-B), E{B-V), E{V-R), E{V-I), 
E{V-H) and E{V~K) against E{V - J) in Fig 8. For nor- 
malization we have used the colour excess EiV — J) instead 
of E{B—V), because E{V—J) does not depend on properties 
like the chemical composition, shape, and structure, degree 
of alignment of the interstellar matter (Sagar & Qian 1990). 
Also, they are better measure of the total amount of inter- 
stellar extinction because of their larger values compared to 
either near-IR or optical colour excesses. A careful selection 
is needed in the case of colour excesses, because the young 
clusters are embedded in emission nebulosity and also con- 
tains young stellar objects. In such environment, the blueing 
effect, ultra-violet excess, circumstellar dust and gas shells 
etc, may be present in and around the cluster stars. So, it is 
fruitful to use V band rather than U or B, because it is least 
affected in such cases. Similarly, in the near-IR, J band is 
preferred, with a view to minimizing the contributions from 
the possible presence of circumstellar material etc, around 
young stellar objects and also to choose a photometric band 
which most closely represents the emission from the stellar 
photosphere. We have therefore estimated the colour excess 
ratios relative to E{V — J). 

In the Fig 8, solid line represents the least square linear 
fits to the data points. The values of correlation coefficient 
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Table 10. Spatial variation of E{B — V) across the cluster Tr 1 and Be 11. The mean values of E{B — V) 
with their standard deviation in mag in 1' X 1' areas are indicated in the appropriate boxes, with the number 
of stars used for this purpose given in brackets. The Aa and AS values are in arcmin relative to the cluster 
center given in Table 2. 



A(5 i 


-3 to -2 


-2 to -1 


Tr 1 

-1 to 


to 1 


1 to 2 


-3 to -2 




0.60 












(1) 








-2 to -1 


- 


0.60 


- 


0.73±0.21 


- 






(1) 




(3) 




-1 to 


0.66 


0.61 


0.60±0.09 


0.58 


0.53 




(1) 


(1) 


(5) 


(1) 


(1) 


to 1 


- 


0.67±0.07 


0.65±0.04 


0.85±0.20 


0.67±0.07 






(2) 


(3) 


(2) 


(3) 


1 to 2 


- 


0.68 


- 


- 


0.64 






(1) 






Cll 


2 to 3 


- 


0.79 


0.71 


- 








(1) 


(1) 












Be 11 






-2 to -1 






1.01±0.02 


l.OOitO.Ol 










(2) 


(2) 




-1 to 






0.95 


0.95±0.05 


1.03 








(1) 


(2) 


(1) 


to 1 




1.05±0.20 


0.98±D.04 


0.94±0.05 


0.99±0.01 






(2) 


(5) 


(3) 


(2) 


1 to 2 












2 to 3 








0.94 





(1) 



(r) and fit indicate that the data points are well represented 
by linear relation. The slopes of these straight lines represent 
reddening directions in the form of colour excess ratios are 
given in Table 11. For comparison, the colour excess ratios 
given by Cardelli et al. (1989) are also listed in the Table 
11. The present reddening directions agree well with those 
given for normal interstellar extinction law. 

In the Fig 9, we have plotted ratios E{V - K)/ E{B -V) 
and E{V - K)/E{V - J) against E{V-K). The horizontal 
line in the fig represents the value of the ratio for the normal 
interstellar extinction law. In the case of normal extinction 
law, the ratio E{V - K)/E{B - V) and E{V - K)/E{V ~ J) 
remains the same for all the values of E{V — K) (Sagar & 
Qian 1990). Least square linear fits to the data points gives 



E{V - K)/E{B -V) = 0.10(±0.03)£;(V^ - K) + 

2.13(±0.31) r=0.29 
E{V - K)/E{V - J) = 0.03(±0.02)i;Cl/ - K) + 

1.18(±0.05) r=0.44 



E(V-K) 



The values of r are < 0.5. This indicates that the relations 
are not statistically significant. This may therefore imply 
absence of anomalous interstellar extinction law toward the 
cluster Be 11. 

Whittet & Breda (1980) suggested that, in the absence 
of complete data at long wavelengths, the approximation 
R = 1.1E{V - K)/E{B - V) generally used to deduce R 
is relatively insensitive to the reddening law adopted. We 
have therefore used this relation to evaluate R. The average 
value of R = 2.60 ±0.34 (sd), which is not too different 
from the value 3.1 for normal extinction law. In the light of 



Figure 9. The colour excess ratios E{V — K)/E{V — J) (top 
panel) and EiV — K)/E{B — V) (bottom panel) as a function 
of E{V — K) for Be 11. The horizontal lines are drawn for the 
normal values of the colour excess ratios. 



above analysis, we conclude that interstellar extinction law 
is normal towards Be 11 in agreement with our earlier result. 
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4.4 Distance to the clusters 

The distances of the clusters are derived by a ZAMS fitting 
procedure. Wo have plotted intrinsic CM diagrams for Tr 1 
and Be 11 in Fig 10. In order to reduce field stax contaminar- 
tion, we have used only those probable cluster member stars 
which are within the cluster radius from the cluster center. 
For plotting these diagrams, we have converted apparent V 
magnitude and {U - B), {B - V), (V - R) and (V - I) 
colours into intrinsic one using the available individual val- 
ues of E{B — V) and following relations for E(U — B) (cf. 
Kamp 1974; Sagar & .Joshi 1979), A„ and E{V - 1) (Walker 
1987) and E{V - R) (Alcala' et al. 1988). 

E{U - B) = [X + Q.QhE{B - V)]E{B - V) 

where X = 0.62 - 0.3(B - V)o for {B - V)o < -0.09 
and X = 0.66 + 0.08(B - V)„ for(B ^ V)o > -0.09 

A^ = [3.06 + 0.25(B - V)o + 0.05E{B - V)]E{B - V); 
and E{V -R) = [El + E2E(B - V)]E{B - V) 

where El = 0.6316 + 0.0713(5 - V)o 
and E2 = 0.0362 + 0.0078(5 - V)o; 

E(V - I) = 1.25[1 + 0.06(B - V)o + 0.014.E{B - 
V)]E{B - V) 

For fainter stars, the average E(B — V) value have been 
used for both the clusters since individual values are not 
known. 

In Vb, (U — B)o and Vo, {B — V)o diagrams, we fitted 
the ZAMS given by Schmidt-Kaler (1982) while the ZAMS 
given by Walker (1985) was fitted in Vo, {V — I)q diagram. 
For Vb, {V — R)o diagram, we have calculated {V — R)o using 
its relation with (B — V)o given by Caldwell et al. (1993). 
The visual fit of the ZAMS to the bluest envelope of the CM 
diagrams gives the mean values of (m — M)o as 12.1±0.2 and 
11.7±0.2 mag for Tr 1 and Be 11 respectively. The distances 
to the clusters should be considered reliable because they 
have been derived by fitting the ZAMS over a wide range 
of the cluster MS. The distance modulus determined above 
yields a distance of 2.6±0.10 Kpc to Tr 1 and of 2.2±0.10 
Kpc to Be 11. For Tr 1, our values of distance modulus is 
in good agreement with the value 12.10 given by Phelps & 
Janes (1994 ). However, it is slightly higher than the value 
11.64 mag given by Joshi & Sagar (1977). The present deter- 
mination of the distance modulus to the cluster Be 11 agrees 
very well with the value of 11.7 mag given by Jackson et al. 
(1980) . 



4.5 Gaps in MS 

The intrinsic CM diagrams of both the clusters Tr 1 and Be 
11 exhibit gaps at different points in the MS (see Fig 10). 
There seems to be a gap between 10.7 and 11.6 intrinsic 
V mag in Tr 1. Another feature seen in this cluster MS is 
the deficiency of stars between 14.0 and 14.8 mag. In the 
case of Be 11, we noticed a gap between 12.4 and 12.9 in- 
trinsic V mag. The reality of the gaps in MS is tested by 
the method adopted by Hawarden (1971). The probability 



of finding such gaps to be accidental is 0.3% for the gap be- 
tween 10.7 and 11.6 in Tr 1 and 0.02% for the gap between 
12.4 and 12.9 in Be 11. The accidental probability is 12% for 
the fainter gap in Tr 1. The very low values of probabilities 
indicate that the observed gaps are real. However, cause of 
such gaps are not well understood. 

4.6 Ages of the clusters 

The ages of the clusters namely Tr 1 and Be 11 have been 

determined by fitting the theoretical stellar evolutionary 
isochrones given by Schaller et al. (1992) in the correspond- 
ing CM diagrams (Fig 10). The isochrones are for Pop I 
stars (X = 0.70, Y = 0.28, Z = 0.02) and include the effects 
of mass loss and convective core overshooting in the model. 
The isoclirone fitting to the main sequence and brighter stars 
indicates that ages of the clusters Tr 1 and Be 11 are 40±10 
and llOzblO Myr respectively. The present age estimate for 
Tr 1 and Be 11 arc in good agreement with the values given 
in Mermilliod (1995) (sec Table 1). For Tr 1, Joshi & Sagar 
(1977) have given the age 3x10^^ yr, which is also in good 
agreement. 

To determine the age and distance of the cluster Be 
11 with the combination of optical and Near-IR data, we 
plot V vs {V - K) and K ys {J - K) in Fig 11 . We 
have overplotted the theoretical isochrones of log(agc) = 
8.05 given by Schaller et al. (1992). The apparent distance 
moduli {m — M)v,(v-K) and {m — MJuxj-K) turn out to be 
14.6±0.3 and 12.0±0.3 mag respectively. By using the red- 
dening estimated in the previous section we derive a distance 
of 2.1±0.3 Kpc for Be 11. Both age and distance determi- 
nation for Be 11 are thus in excellent agreement with our 
earlier estimates. 

4.7 Luminosity and Mass function of the clusters 

The luminosity function (LF), denotes the relative number 
of stars in unit magnitude range. The correction of non- 
member stars is very much important in the construction 
of LF's for star clusters. Two colours, such as B and V or 

V and I, are required for the non-member identification. It 
is therefore required to construct the LF cither from a V, 
{B — V) diagram or from V, {V — I) diagram or from a 
similar diagram instead of from a single B, V or I or any 
other passband. We preferred the V, {V — I) over the other 
diagrams as it is deepest. The main disadvantage of using 
two passbarids for the construction of the LF is that both 
passbands introduce incompleteness, whose determination is 
a difficult process, as described below. 

4-7.1 Determination of photometric completeness 

The method consists of insertion of randomly selected artifi- 
cial stars with known magnitude and position in the original 

V frame. For the I band image the inserted stars have same 
geometrical positions but differ in I brightness according 
to mean {V — I) colour of the MS stars. Only 15% of the 
number of actually detected stars arc inserted at one time, 
so that the crowding characteristics of the original data re- 
mains almost unchanged. The luminosity distribution of the 
artificial stars has been chosen in such a way that more stars 
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Table 11. A comparison of extinction law in the direction of Be 11 with normal extinction law given by 
CardeUi et al. (1989). 



Source 


E(U-B) 
E(V-J) 


E{B-V) 
E(V-J) 


E{V-R) 

E(y-J) 


E(V-I) 
E{V-J) 


E{V-H) 
E{V-J) 


E(V-K) 
E(V-J) 


E(J^K) 
E{V-K) 


Cardelli et al. 
Colour excess ratio 


0.32 
0.42±0.13 


0.43 
0.53±0.03 


0.27 
0.29±0.03 


0.56 
0.56±0.04 


1.13 
1.31 ±0.07 


1.21 
1.38±0.07 


0.19 
0.18±0.30 




-10 1 12 10 1 

{U-B)„ (B-V)„ (V-R)o (V-I)o 



Figure 10. The Vo, {U - B)o; Vb, {B - V)o; Vo, {V - R)o and Vb, (V - I)o diagrams for stars of the Tr 1 and Be 11. The continuous 
curves are the ZAMS fitted to the MS. The dotted curves are the isochrones for Pop I stars of log(age)=7.6 for Tr 1 and 8.05 for Be 11 
fitted to the brighter cluster members. 



are inserted into the fainter magnitude bins. The frames are 
re-processed using the same procedure used for the original 
frames. The ratio of recovered to inserted stars in the dif- 
ferent magnitude bins gives directly the completeness factor 
(CF) for that region. Table 12 lists the CF in both cluster 
and field regions of the objects under study. 

Several authors have discussed the method for deter- 
mining the CF (cf. Stetson 1987; Mateo 1988; Sagar & 



Richtler 1991; Banks et al. 1995). Sagar & Richtler (1991) 
have taken the minimum value of the completeness factors 
of the pair to correct the star counts. They argue that the 
two frames are not independent and that the multiplicative 
assumption of Mateo (1988) could not be justified. Banks et 
al. (1995) tested the ability of the techniques given by Mateo 
(1988) and Sagar & Richtler (1991) on the basis of numerical 
simulations. They conclude that product method of Mateo 
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Figure 12. The plot shows the mass functions derived using 
Schaller et al. (1992) isochrones. 



Figure 11. The K, (J - K) and V, {V - K) diagrams of the 
sample stars in the cluster Be 11. 



Table 12. Variation of completeness factor (CF) with the MS 
brightness in both cluster and field regions. 



Tr 1 Be 11 



V mag 


range 


Cluster 


Field 


Cluster 


Field 


13 - 


14 


0.96 


0.99 


0.99 


0.99 


14 - 


15 


0.96 


0.99 


0.99 


0.99 


15 - 


16 


0.95 


0.99 


0.99 


0.99 


16 - 


17 


0.94 


0.98 


0.98 


0.99 


17 - 


18 


0.94 


0.98 


0.97 


0.99 


18 - 


19 


0.90 


0.98 


0.96 


0.98 


19 - 


20 


0.80 


0.98 


0.93 


0.98 



(1988) increasingly over-estimates the incompleteness cor- 
rection as the magnitude is increased, and the method sug- 
gested by Sagar & Richtler (1991) recovered the actual LF 
better with a mean error of 3% upto CF > 0.5. We therefore 
used the procedure of Sagar & Richtler (1991) in the present 
work. 

4-7-2 Determination of Mass Function 

To derive the true LF of the cluster we have to remove the 
field star contamination. For this, we use photometric crite- 
rion. First, we defined a blue and red envelope for the MS 
on V , (y — I) CM diagram. Same envelope were drawn on 
the corresponding field region. Star counts as a function of 
luminosity were made in both (cluster and field) regions. In 
this way we can estimate the number of field stars present in 
various magnitude bins of the cluster region. The observed 
LFs of the cluster and field regions were also corrected for 
data incompleteness as well as for differences in area. True 
LF for the cluster was obtained by subtracting the observed 
LF of field region from the observed LF of cluster region. 

The mass function (MF), which denote the relative 
number of stars in unit range of mass centered on mass M. 



Table 13. The slope of the mass function derived from LF along- 
with relaxation time T^. 



cluster 


Mass range 


Mass Function slope 


log Te 




Mo 


{x) 




Tr 1 


0.9 - 5.1 


1.50±0.40 


7.2 


Be 11 


1.0 - 4.5 


1.22±0.24 


7.1 



The MF slope has been derived from the mass distribution 
5(M). If dN represents the number of stars in a mass bin dAI 
with central mass M, then the value of slope x is determine 
from the linear relation 

log|^ = -(l+a;)xlog(Af)+constant 

using the least-squares solution. The Salpeter (1955) value 
for the slope of MF is a; = 1.35. 

To derive the MF from LF, we need theoretical evolu- 
tionary tracks and accurate knowledge of cluster parame- 
ter like reddening , distance, age etc. Theoretical models by 
Schaller et al. (1992) were used to convert the observed LF 
to the MF. Fig 12 represents the plot of MFs of Tr 1 and 
Be 11. The value of the MF slope along with the mass range 
and error are given in Table 13, where the quoted errors are 
errors resulting from the linear least square fit to the data 
points. For the cluster Tr 1, the value of x is in agreement 
with the value given by Phelps & Janes (1993). The values 
of X for both clusters are in agreement with the Salpeter 
value. 

The dip is seen in the MF of both the clusters. In 
the cluster Tr 1 one dip is present in the mass range 
1.5< M/Mq <1.9. This dip is due to the deficiency of stars 
in the intrinsic V magnitude between 14.0 and 14.8 (see Sec. 
4.5). In this cluster an apparent dip in LF is also found by 
Phelps & Janes (1993) at 1.4 Mq- One dip is also seen in the 
MF of Be 11 in the mass range 2.0< M/Mq <3.0. This dip 
corresponds to the gap noticed in the intrinsic V magnitude 
between 12.4 and 12.9 (see Sec. 4.5). 
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5 MASS SEGREGATION 



To study the eflect of mass segregation for the clusters under 
study, we plot in Fig. 13 cumulative radial stellar distribu- 
tion of stars for different masses. A careful inspection of Fig 
13 shows that both the clusters have mass segregation effect. 

To perform the K-S test among these distribution to 
see whether they belong to the same distribution or not, 
we have divided in three mass range 5.0< M0<2.5, 2.5< 
Mq<1.0 and M0<1.O and 4.5< Mo<2.5, 2.5< M0<1.5 and 
M0<1.5 mag for Tr 1 and Be 11 respectively. The K-S test 
indicates that mass segregation has occurred at confidence 
level of 99% for Tr 1 and 80% for Be 11. One would like to 
know whether existing mass segregation is due to dynamical 
evolution or imprint of star formation process. 

Dynamical evolution is one of the possible cause for 
mass segregation. At the time of formation, the cluster may 
have a uniform spatial stellar mass distribution, which may 
be modified due to dynamical evolution of the cluster mem- 
bers. Because of dynamical relaxation, low mass stars in a 
cluster may possess largest random velocities, consequently 
these will try to occupy a large volume than the high mass 
stars (cf. Mathieu & Latham 1986, McNamara & Sekiguchi 
1986, Mathieu 1985). Thus mass segregation develops in the 
time scale required to exchange energy between stars of dif- 
ferent mass by scattering. The dynamical relaxation time, 
T_B is the time in which the individual stars exchange ener- 
gies and their velocity distribution approaches a maxwellian 
equilibrium. It is given by 



Te 



8.9 X W^N^/'^rI^^ 



< m >i/2 log(OAN) 

where A'' is the number of cluster members, is the 
radius containing half of the cluster mass and < m > is the 
average mass of the cluster stars (cf. Spitzer & Hart 1971). 
The number of probable MS stars is estimated using the 
CMDs of the clusters after subtracting the contribution due 
to field stars and applying the necessary corrections for the 
data incompleteness. Due to our inability to estimate the 
from the present data, we assume that the R^ is equal to half 
of the cluster radius derived by us. The angular values are 
converted to linear values using the cluster distances which 
are derived here. Inclusion of cluster members fainter than 
the limiting V magnitude will decrease the value of < m > 
and increase the value of A''. This will result in higher values 
of Te ■ Hence the Te values obtained here may be considered 
as the lower limit. 

A comparison of cluster age with its relaxation time 
indicates that the relaxation time is smaller than the age of 
the clusters. Thus we can conclude that the clusters under 
study are dynamically relaxed. It may be due to the result of 
dynamical evolution or imprint of star formation processes 
or both. 



6 CONCLUSIONS 

Physical parameters of the young open star cluster Tr 1 and 
Be 11 have been derived based on UBVRI CCD photome- 
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Figure 13. Cumulative radial distribution of stars in different 
mass ranges for Tr 1 and Be 11. 



try. Our study leads to the following conclusions: 

(i) The radial density profiles of the clusters indicate the 
existence of clustering. The angular diameter of the cluster 
Tr 1 and Be 11 are 6'.0 and 4'. 8 respectively. The corre- 
sponding linear sizes are 4.6 and 3.0 pc respectively. 

(ii) Variable reddening is present within both the clus- 
ters Tr 1 and Be 11 with the mean value of E{B — 1/) = 
0.60±0.05 mag and 0.95±0.05 mag respectively. The law 
of interstellar extinction is normal in the direction of both 
the clusters. Combining 2MASS data with optical data, we 
have studied the extinction law in the direction of Be 11. 
Colour-colour diagram yields the colour excess i7( J — K) = 
0.40±0.20 mag and E{y-K) = 2.20±0.20 mag respectively. 

(iii) The distance values are 2.6±0.10 and 2.2±0.10 Kpc 
for the cluster Tr 1 and Be 11 respectively. The fitting of 
Schaller et al. (1992) isochrones to the intrinsic CM dia- 
grams indicate an age of 40±10 and 110±10 Myrs for for Tr 
1 and Be 11 respectively. 

(iv) The luminosity function of clusters was constructed 
by subtracting field star contamination determined from 
neighbour field's. The luminosity function was transformed 
into the mass function using the present cluster parameters 
and the theoretical model given by Schaller et al. (1992). 
This gives the mass function slope x — 1.50 ± 0.40 and 
1.22±Q.24 for Tr 1 and Be 11 respectively. The values of 
X for Tr 1 and Be 11 thus agrees with the Salpeter (1955) 
value. 

(v) Mass segregation, in the sense that massive stars 
tend to lie near the cluster center, is observed in both the 
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cluster Tr 1 and Be 11. The dynamical relaxation time in- 
dicate that both the clusters are dynamically relaxed and 
mass segregation may have occurred due to dynamical evo- 
lution, or imprint of star formation or both. 

(vi) Present data are unable to unambiguously identify 
the members of the clusters. For this, precise kinematical 
observations are required. 
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